Charging effects in the inductively shunted Josephson junction 
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The choice of impedance used to shunt a Josephson junction determines if the charge transferred through 
the circuit is quantized: a capacitive shunt renders the charge discrete, whereas an inductive shunt leads to 
continuous charge. This discrepancy leads to a paradox in the limit of large inductances L. We show that while 
the energy spectra of the capacitively and inductively shunted junction are vastly different, their high-frequency 
responses become identical for large L. Inductive shunting thus opens the possibility to observe charging effects 
unimpeded by charge noise. 

PACS numbers: 85.25.Cp, 74.50.+r, 03.65.-w 



The Josephson junction plays a key role in superconducting 
devices, being the only nonlinear yet dissipationless element 
in a quantum circuit [l]. The simplest examples of nonlinear 
quantum circuits are obtained by shunting a single Josephson 
junction by a purely dispersive impedance Z(ui). Following 
this prescription and using a capacitive shunt, one obtains the 
Cooper pair box (CPB) Jldll, see Fig. da). Alternatively, 
using an inductive shunt, one obtains the single-junction flux 
qubit |5[ and the phase qubit see Fig.[TJb). 

There is an important difference between these two exam- 
ples: in the circuit of Fig. [TJa), one of the terminals of the 
Josephson junction forms a superconducting island coupled 
only capacitively to the rest of the circuit, while the circuit of 
Fig. 02b) does not contain isolated pieces of superconductor. It 
is well-known that charges of two superconductors coupled by 




FIG. 1: (color online) Examples of dispersively shunted Josephson 
junctions, (a) Island-based devices, like the Cooper pair box, are 
obtained by shunting a Josephson junction capacitively and are typi- 
cally operated by coupling to charge, (b) Loop-based devices like the 
one-junction flux or phase qubit use an inductive shunt of the Joseph- 
son junction and are commmonly operated by coupling to flux, (c) 
For large inductances as in the fluxonium device, both charge and 
flux coupling play a role, (d) According to Thevenin's theorem, 
both circuits (a) and (b) represent specific realizations of a Joseph- 
son junction shunted by an impedance Z(uj). Taking the limit of 
large inductance should smoothly connect the realization (c) to (a). 



only a Josephson junction are quantized at the level of eigen- 
values of the corresponding operators. Adding an inductive 
link between the superconductors destroys that quantization. 

Paradoxically, intuition suggests that charge quantization 
should manifest itself as long as the inductance remains suf- 
ficiently large. Specifically, the properties of the inductively 
shunted junction [Fig.QJc)] should approach those of the CPB 
[Fig. QI a )] as the inductance L is increased. Indeed, us- 
ing Thevenin's theorem, the environment of the Josephson 
junction can formally be described by the series impedance 
Zj h = iu)L(\ — up-LCg)' 1 [Fig. [Ttd)], which converges to 
the CPB expression Zcpb = (* W C S ) _1 as L — > oo. 

The central question formulated and answered in this Let- 
ter is whether and in what sense the charge of finite-size su- 
perconductors remains quantized in the presence of a material 
link between them. We find the proper observable quantity 
that allows one to see the manifestation of charging effects 
and its evolution with the inductance of the material link. In 
addition to the fundamental interest, our observation has a 
practical value Q]: excited states of the inductively shunted 
Josephson junction turn out to be robust with respect to charge 
fluctuations while being strongly anharmonic, which is an ad- 
vantage over earlier existing superconducting devices. 

We start with the Hamiltonian describing the inductively 
shunted junction as shown in Fig. [He), 

H = AE c {n~n a f-Ej C os V + ^E L { V + 2^/^ ) 2 . (1) 

Here, the canonically conjugate operators n = q/2e and ip 
correspond to the junction capacitor charge (in units of Cooper 
pairs) and the superconducting phase difference across the 
junction, and they obey the relation [ip, n] = i. The three en- 
ergy scales entering H are the (single-electron) charging en- 
ergy Eq = e 2 /2(Cj + C g ), the Josephson energy Ej, and the 
inductive energy El = {^o/2tt) 2 / L, defined in terms of the 
flux quantum $o = h/2e. The external magnetic flux through 
the loop (generated by the bias current /;,) is denoted by <5, 
and the effect of the gate voltage is expressed in the offset 
charge n g . We note that the spectrum is invariant with respect 
to static offset charge, which can be confirmed by applying a 
gauge transformation \[){<p) = e" lg '- p ^j(ip). 

In spite of the familiarity of the Hamiltonian (Q]), its physics 
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beyond the flux and phase qubit regimes (Ec <C El < Ej) 
has remained largely unexplored, with the notable exception 
of the work by A. Kitaev |8|. We focus on the regime of 
large inductances where El represents the smallest energy 
scale, El <C Ec, Ej. In the first step, we show that large- 
inductance limit does not smoothly transform the energy spec- 
trum of the inductively shunted device into that of the CPB. 
Ultimately, the reason for the subtlety of this limit stems from 
the opposed symmetries of the two systems: while the CPB 
has a strictly periodic potential when written in the phase ba- 
sis, the inductive shunt always breaks this periodicity. Equiv- 
alently, the discrepancy may be formulated in the charge ba- 
sis: due to the presence of an island in the CPB, charge on 
the corresponding node is quantized in terms of Cooper pairs, 
whereas the inductive shunt eliminates the island and renders 
the charge continuous. 

The key to our analysis consists of transforming the Hamil- 
tonian (Q]) into the basis of Bloch waves {|s,p)}, where seN 
is the band index and p G [0,1) the quasimomentum. These 
states diagonalize the CPB part [comprising the terms ~ Ec 
and Ej] of the Hamiltonian (Q~|), Ho\s,p) = e s (p)\s,p), 
where the eigenenergies represent the usual CPB bands. For 
the transformation of the inductive term, we employ the stan- 
dard relation ip = id/dp+Cl jjj, where the Hermitian operator 
f2 causes interband coupling and is defined via 



(ps\n\p's') = S(p-p'} 
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Here, u ps is the Bloch amplitude, ((p\s,p) 
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large Ej /Ec and low-lying bands s, s' we find fl ss > (p) « 
(2E c /Ej) 1 / 4 (y/sS StS > + i + Vs 7 <5 s , S '-i)/27r so that interband 
coupling can be neglected. This simplifies the problem dra- 
matically: the Hamiltonian becomes block-diagonal and sep- 
arates into effective Hamiltonians 
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for the low-lying bands. The spectra, obtained from H ( s > with 
periodic boundary conditions in p, have to be overlaid. We 
note that Eq. ([3]l is structurally identical to the CPB Hamilto- 
nian, and in the following we will draw from the knowledge 
of the CPB in both charge and transmon regimes SOU- 

According to Eq. (01, El determines the kinetic energy 
while the CPB bands e s (p) act as potentials for the induc- 
tively shunted device. The potentials are periodic and can be 
expressed as Fourier series e s (p) — J2'eLo £s < 1 cos(2tt£p), and 
approach simple sinusoids for large Ej / Ec- For small induc- 
tive energies, each low-lying CPB band supports two types of 
states: (i) metaplasmon states with approximate level spac- 
ing 2-k^JEl which are bound states in the CPB poten- 
tial, and are associated with large charge oscillation across the 
Josephson junction, and (ii) persistent-current states, existing 
in the gaps of the "bare" (El = 0) CPB device; these states 
are associated with substantial persistent currents circulating 
(at finite El) through the inductor. As a result, the energy 
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FIG. 2: (color online). Energy spectra for the inductively shunted 
losephson junction at zero flux. The progression (a)-(c) shows en- 
ergy spectra for Ej/Ec = 2.5 and decreasing inductive energy 
El. Energies are given in units of the plasma oscillation frequency 
tiujp = \/8EjEc- For small El, the spectrum exhibits regions with 
significantly different level spacings. (d) Comparison with the CPB 
bands for the same Ej/Ec ratio reveals that regions of small (large) 
level spacings coincide with band (gap) regions in the CPB. 



spectrum of the shunted junction separates into regions with 
small level spacings fitting into the bands of the CPB, and 
large level spacings in the regions of CPB gaps. This distinc- 
tion grows stronger with the decrease of the ratio El /Ec, see 
Fig.HaMc). 

The nature of the persistent-current states can further be 
clarified by rewriting the Hamiltonian (O in the discrete local- 
minimum basis {|m, s)}, which formally parallels the charge 
basis of the CPB, 



H 



(s) _ 
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\m, s) (m + 1, s | + H.c. 



Physically, |m, s) corresponds to a state localized in the m-th 
local minimum of the phase-basis potential, carrying a persis- 
tent current 1^ = (to$o + < E > )/i. From Eq. ©, the energies 
of persistent-current states are found to be 

E (s) 



(5) 
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where m £ Z has sufficiently large modulus so that E, 
maxp e s (p). Thus, to lowest order in the inter-well tunneling, 
the persistent-current states are doubly degenerate at zero flux, 
and each degenerate subspace is spanned by the two counter- 
propagating states with currents ~ ±m$o / L around the 
superconducting loop. At higher order in the inter-well tun- 
neling, this degeneracy is lifted and time-reversal symmetry, 
which requires 1 = for vanishing magnetic flux, is restored. 
The resulting avoided crossings are extremely small. Specif- 
ically, for a purely sinusoidal s = band, i.e. Eqj> = for 
I > 1, the splitting between \m,s = 0) and |— m, s = 0) is 
generated by tunneling through the 2m potential barriers sep- 
arating the two states, and scales as 
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FIG. 3: (color). Metaplasmon and persistent-current states for large 
inductances. Panel (a) shows the potential in i^-space (cyan), the 
corresponding bent CPB bands (orange), and examples of wavefunc- 
tions at zero flux, from bottom to top: black-ground-state wavefunc- 
tion (lowest s=0 metaplasmon state), blue-s=0 persistent-current 
states, green-lowest s=l metaplasmon state, and yellow-s=2 meta- 
plasmon state. Localization in <^3-space differs characteristically, 
with metaplasmon states being centered at <p = and delocal- 
ized, while persistent-current states are symmetric and antisymmet- 
ric superpositions of wavefunctions localized towards the edges of 
the parabolically deformed bands, (b) Corresponding CPB spec- 
trum for comparison. Using the same color coding, (c) and (d) show 
the spectrum overlayed by the CPB bands, and the flux dependence 
of energy levels, respectively. Metaplasmon states (independent of 
flux) and persistent-current states (sensitive to flux) are easily dis- 
tinguished. For the parameters chosen here, Ej/Ec = 2.5 and 
El/Ec = 10 -3 , interband coupling becomes significant at higher 
energies (levels in gray), leading to avoided crossings between s = 1 
persistent-current states (in red) and s = 2 metaplasmon states. 



Metaplasmon and persistent-current states also differ in 
their localization and respective flux dependence. As evident 
from Eq. and illustrated in Fig. a), metaplasmon states 
are typically delocalized across several wells around ip = 0, 
whereas persistent-current tend to localize in </>space in the 
region of the corresponding parabolically deformed band. For 
the bound metaplasmon states, the flux dependence is ex- 
ponentially suppressed; specifically, variations in the meta- 
plasmon states of band s scale as ~ exp[ 



By contrast, Eq. (0 predicts that persistent-current states are 
strongly sensitive to magnetic flux. For the lowest bands, 
these predictions are confirmed by results from numerical di- 
agonalization, see Fig. [3jd). Persistent-current states exist 
above the lowest CPB band and overlap with s > metaplas- 
mon states. At higher energies, avoided crossings between 
metaplasmon and persistent-current states become more sig- 
nificant, as interband coupling increases for higher bands. 

The previous discussion underlines the differences between 
the spectra of the CPB and the inductively shunted Josephson 
junction: the presence of the inductor introduces new levels 
located in the band gaps of the CPB, and these levels do not 
disappear in the limit of small E^. However, we will now 
argue that in the limit of large inductance, the ac properties 
of the shunted device approach those of the CPB. To be spe- 
cific, we consider the realistic ifioll coupling of the junction 
charge n to an ac voltage potential, which modifies the off- 
set charge term in Eq. (1), n g — > n g + C g V ims (a + a^)/2e. 
Here, a, annilate or create a photon in the microwave field 
mode, Vrms denotes its root-mean square voltage. Transitions 
are then induced by tuning the field into resonance, and the 
transition strength depends directly on the charge matrix el- 
ements (f\n\i), where i, f denote the initial and final states 
of the transition. At low energies where interband coupling is 
negligible, these states can be identified as eigenstates of Eq. 
(01, and we will write, e.g., \i) = |s„), thus denoting the v-th 
eigenstate belonging to band s. We find 
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{s' vl \n\s v ) =i(Ej/2Ec) 1/4 (V^S SiS , + i - VdS 

,1/2 

dpXs'v'(p)Xsv(p)- (7) 

1-1/2 

The last integral involves the quasimomentum representation 
Xsu{p) = (sp\s v ) of eigenstates, and is amenable to eval- 
uation within the WKB approximation. In complete anal- 
ogy to the Franck-Condon principle, the matrix elements 
peak when the transition between the corresponding classi- 
cal turning points occurs vertically, see Fig. [4] For the meta- 
plasmon transitions from s = to 1 we find that below 
this optimum the matrix elements decay exponentially with 
~ exp[— y/2 |eo,i| E~l(w/ |£i,i|) 3 ' 2 /27t], whereas for ener- 
gies above, the decay is oscillatory and follows an overall 
power-law decay ~ [£ , |/(|e 11 e 01 | rj)] 1 / 4 , where r\ is the en- 
ergy deviation from the optimum. As a result, in the limit of 
large inductances the matrix elements select the CPB-allowed 
transitions, thus producing pronounced charging effects. 

We now discuss the charge-noise sensitivity of the induc- 
tively shunted junction, and show that the charging effects re- 
vealed in the ac response remain largely unspoiled by 1/f 
charge noise. The origin of the reduced sensitivity to charge 
noise is identical to that of the the single-junction flux qubit 
and the flux-biased phase qubit 111 111 . Closing the Josephson 
junction by an inductive load allows one to evade the most se- 
vere problems caused by 1 // charge noise. Specifically, the 
vanishing dc impedance of the inductive element renders the 
energy spectrum independent of any constant offset charge, 
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FIG. 4: (color online). Matrix elements relevant for excitation when coupling the junction charge to an ac voltage. The magnitude of the 
matrix element for a transition to final state \j) is shown as a function of the final state energy Ej, in (a) for El/ Ec = 1CP 3 starting in the 
ground state, and in (b) for El/ Eq = 10~ 4 for several initial states. The prominent gap in the matrix elements (right panels) demonstrates the 
strong suppression of transitions from metaplasmon states to persistent-current states. For transitions among metaplasmon states, the matrix 
elements reach a distinct maximum for the transition occuring in the corresponding CPB without the inductive shunt (vertical dashed lines in 
left panels). Lowering El makes the peak narrower, and eventually leads to a selection rule excluding all but the CPB-allowed transitions. 



and transforms the 1// charge noise into a relatively benign 
"/-noise". 

To confirm this, consider the effects of a fluctuating off- 
set charge n g (t) in the Hamiltonian (Q~|i. Again, we apply a 
gauge transformation V>(<^) = e" lsip ip(ip). Taking into ac- 
count the time dependence of n g , this results in a Hamilto- 
nian H, obtained from H by eliminating the offset charge 
n g from the charging term and appending a term —Hh g (p. 
Thus, the energy spectrum only depends on the time deriva- 
tive of the offset charge, in a fashion similar to the depen- 
dence on flux see Eq. (Q3 lfl2l] . The relevant noise spec- 
trum for fig is Sn g (uj) — Lo 2 S ng (uj), such that Gaussian 
charge noise with a 1// spectrum, S„ (ui) = A/ \ui\, is trans- 
formed into benign "/-noise". Assuming an ultraviolet cut- 
off of the 1// spectrum at frequency l/r u , the residual ef- 
fect of n g noise away from flux sweet spots and for long 
times t ^S> t u is controlled by the behavior of the integral 
g(t) = (duiij j dhg) 2 fijp 1 did S„ g (ui), such that off-diagonal 
elements of the density matrix decay as pij ~ exp[— g(t)]. For 
large t, the integral g(t) is logarithmically divergent, and thus 
yields a slow power-law decay with exponent proportional to 
{duiij / '9$) 2 . Remarkably, the flux sweet spots thus become 
first-order insensitive to charge noise, and metaplasmon states 
with their suppressed flux dependence are expected to be es- 
pecially well-protected. The slow power-law decay should be 
contrasted with the rapid loss of coherences in the absence of 
an inductive shunt, which follows an exponential decay. 

To summarize, the physics of the inductively shunted 
Josephson junction in the large-i limit is distinct from the 
ones accessed by flux and phase qubits. Two different types 
of states, metaplasmon and persistent current states, with dis- 



tinct level spacings and magnetic flux dependence dominate 
the spectrum at low energies. While the spectrum of the junc- 
tion with large inductive shunt also differs from that of the 
CPB, we have demonstrated that the ac response due to charge 
coupling approaches the well-known response of the CPB, 
thus resolving the dichotomy of island vs. loop based devices 
in the L — > oo limit. These findings have been successfully 
employed in the analysis of a recent experiment 10], and will 
be of future interest in exploring the device's applicability to 
quantum information processing and observation of Bloch os- 
cillations. 
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